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We study numerically the dynamics of the Rabi Hamiltonian, describing the interaction of a single 
cavity mode and a two-level atom without the rotating wave approximation, subjected to damping 
and dephasing reservoirs included via usual Lindblad superoperators in the master equation. We 
show that the combination of the antirotating term and the atomic dephasing leads to linear asymp- 
totic photons generation from vacuum. We reveal the origins of the phenomenon and estimate its 
importance in realistic situations. 
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I. INTRODUCTION 

A fundamental task in Physics is the description of the 
matter-light interaction. The most simple model to deal 
with is the Rabi model [T], which describes the interac- 
tion of a two-level atom with a single mode of the quan- 
tized electromagnetic (EM) field. The Rabi Hamiltonian 
(RH) reads (fi = 1) 

H = uata + -^ a z + 9 (o+ + cr~) (« + a*) , (1) 

where u> and ujq are the field and atomic transition fre- 
quencies, respectively, and g is the coupling constant 
(vacuum Rabi frequency), a (a?) is the annihilation (cre- 
ation) operator of the EM field, a z = |e) (e| — \g) (g\ 
and cr+ = |e) (g\ (<r_ = 0+) are atomic operators, 
with \g) and |e) denoting the ground and excited atomic 
states, respectively. Although largely studied over the 
last decades, up to now its exact analytical solution is 
lacking and only numerical [2j [3j SJ [5] and approximate 
analytical solutions are available [5J |7J |5] , despite the con- 
jecture by Reik and Doucha [5] that an exact solution of 
RH in terms of known functions is possible. The most 
used analytical approach to RH is to make the rotate 
wave approximation (RWA), where the antirotating term 
g (a'cr+ + aer_) is neglected , since in the weak coupling 
regime g/oj <C 1, small detuning |A| <§; oj (A = ujq — lo), 
and weak field amplitude its contribution to the evolu- 
tion of the system is quite small [TUHH]- In this limit the 
RH is known as Jaynes-Cummings Hamiltonian (JCH) 
[T2l Il3] and can be integrated exactly. 

For having an exact solution, the JCH has been largely 
employed in Quantum Optics, in particular in cavity 
quantum electrodynamics (QED) [HUTS], where the vast 
majority of the experiments satisfy the required parame- 
ters regime [1IJJQT]. JCH revealed interesting phenomena 
related to the quantum nature of the light, encompassing 
the granular nature of the electromagnetic field, revealed 
through collapse and revivals of the atomic inversion [16] , 
Rabi oscillations [T7J, squeezing [13], non-classical states, 



such as Schrodinger cat state-like and Fock states 
|20j . and the entanglement between atom-atom or atom- 
field [21] . The manipulation of atom-field interaction has 
been employed in the implementation of quantum logic 
gates in trapped ions [55] and in cavity QED [23], as 
well as the atomic teleportation process [24], which have 
contributed for a fast development of the quantum infor- 
mation science [2~5] . Moreover, over the past few years 
the JCH was experimentally investigated in solid state 
cavity QED systems in the strong coupling regime using 
superconducting artificial two-level atoms coupled to mi- 
crowave waveguide resonators ^jo\ [571 H3J (the so called 
circuit QED [29 ) and quantum dots coupled to photonic 
crystals microcavities j30j[3T]. In fact, in circuit QED the 
JCH is the basic theoretical tool for describing quantum 
logic gates and read-out protocols [35J [33J . 

The antirotating term neglected under RWA is usually 
wrongly interpreted as being non-conservative |10l 1111 
I3~4"] , once it could allow for a violation of the energy con- 
servation. In fact, this term does not conserve the total 
number of quanta of the system, defined by the operator 
N = a^a + a z + I (where I stands for idendity opera- 
tor). However, it does not change the total energy of 
the system, as can be easily seen through the Heisenberg 
equation of motion for the total Hamiltonian operator 
dH/dt = i [H, H] + dH/dt = 0, once H in Eq. is time 
independent. Therefore, the total energy of the system 
(H) is conserved and no violation of physical laws occur. 
Besides, recent works questioned the validity of the RWA 
[3"51 1551 157] and proposed alternative analytical approx- 
imate methods [3 [S] . Moreover, it was shown that the 
antirotating term is responsible for several novel quantum 
mechanical phenomena, such as quantum irreversibility 
and chaos 38, 39J, quantum phase transitions |40j . imple- 
mentation of Landau-Zener transitions of a qubit in cir- 
cuit QED architecture [HI [42] , generation of atom-cavity 
entanglement [43, 44J, and simulation of the dynamical 
Casimir effect (DCE [15]) in semiconducting microcavi- 
ties HaiHl HE] or circuit QED @3]. 

In this work we study numerically the dynamics of the 
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RH subjected to dissipative effects acting on both the 
atom and cavity mode. These undesirable effects are 
taken into account through the master equation approach 
[10j . where the Lindblad superoperators are built as usual 
[49) . Without any formal prove, we simply assume that 
the evolution of the density operator of the system p(t) 
is described by the master equation 

^ = -i[H,p]+C(p), (2) 

where H is the RH Q and the Lindblad operator C (p) 
is given by 

C(p)=L a (p)+£ d (p)+L f (p), (3) 
with the standard definitions 

£a(p) = 2 («t + 1) (2<7_p<7 + - a + (7-p- pa + a-) 

+ 2 Ut ( 2cr +P cr - - O-V+P - pcr_cr+) , (4) 
£f(p) = 2 ( n t + 1) (2apa^ — a) ap — pa) a) 

+ — n t {2a) pa — aa) p — paa^) , (5) 
£d (p) = Jph {°zPVz - p) ■ (6) 

The superoperators C a (p) and C / (p) describe the ther- 
mal reservoirs effects (with mean photon number n t ) on 
the atom and field, respectively, where 7 (k) is the atom 
(cavity) relaxation rate. Another source of decoherence 
of the atom is the phase damping reservoir, represented 
by Ld (p), where 7^ is the dephasing rate. By focusing 
our attention on the asymptotic photons creation from 
vacuum driven by the combination of the antirotating 
term g (a) <r + + aer_) in (JTj) and the atomic phase reser- 
voir, we show that even in situations where the atom 
and field are initially prepared in their individual ground 
states, i.e., |</>) = \g, 0), where |0) is the ground state of 
the EM field, there is an asymptotic photons generation. 
This process depends on the intensity of the atom-field 
coupling g, the detuning A, and it is considerably am- 
plified when the atomic phase reservoir is predominant 
over the other dissipative channels, such as atomic and 
field damping due to the thermal reservoirs. The essence 
of the photons creation mechanism presented here relies 
on the existence of the antirotating term in the RH and 
the limitation imposed by the quantum vacuum, namely, 
a |0) = 0. The role played by the atomic phase reservoir 
is just to amplify the photons creation process through 
atomic decoherence. As far as we know, this phenomenon 
has not been described in the literature. 

This work is organized as follows. In Sec. H we study 
in details the process of photons generation due to the 
atomic dephasing and explain its origin. In Sec. HI we 
give an alternative physical explanation of photons gen- 
erations by considering that the atom's dephasing is due 
to random shifts of the atomic transition frequency, as 
usually occurs in solid state systems (e.g. circuit QED). 



In Sec. IV we study the influence of damping process on 
the photons creation through decoherence and estimate 
the net effect in realistic situations. Finally, the section 
V contains discussion of results and concluding remarks. 



II. PURE DEPHASING 

For convenience, from now on we set the cavity fre- 
quency to u> = 1 . First we integrated numerically the Eq. 
([2]) considering just the dephasing reservoir (7 = n = 0). 
In Figs, f and 2 we set the parameters A = and g = 0.1. 
These and other values of the parameters were chosen 
in order to optimize the numerical calculations, and are 
not related to the experimental data; nevertheless, we 
checked out that qualitatively the behavior described be- 
low also holds for realistic parameters (see Fig. 5b be- 
low). The state |ff,0) is the ground state of the JCH, 
so it is not coupled to other states under the JCH dy- 
namics. On the other hand, the antirotating term in the 
RH does induce transitions to other states and, since the 
atom and field states are limited from below by \g) and 
|0), respectively, these transitions may only increase (n). 
In Fig. fa we plot the mean photon number (n) without 
dephasing (7^ — 0, line f ) as function of the dimension- 
less time r = £t with £ = 0.1 for the initial state \g, 0), 
showing a bound oscillating behavior as described in de- 
tails in [3]. Fig. fa also shows (n) (line 2) and the Mandel 
factor q = ((An) — (n))/ (n) (line 3) for dephasing rate 
Iph = 0.1 : now (n) increases with time, achieving asymp- 
totically a linear behavior, and the generated field state 
demonstrates a super-poissonian behavior, q > 0. 

The photons creation mechanism through atomic de- 
phasing demonstrated in Fig. la is a general phenomenon 
and asymptotically does not depend on the initial state 
of the system. In Fig. lb we plot (n) versus r for 
6 different initial states: \(f>i) — \g,Q), = Iff, ct) 

where \a) is the coherent state and we take \a\ = 0.05, 
10a) = [(Iff) + |e» /n/2] ® |0), |&) = [(| 5 ) + |e» /y/2] ® 
\ a ): 105) = I e j0), an d |</>6) = \ e : a )- The basic differ- 
ence between these initial states is the amount of quanta 
(N k ) = (cb k \N\cb k ): we have (JVi) = 0, (JV a ) = 0.05, 
(N 3 ) = 0.5, (N 4 ) = 0.55, (iV 5 ) = f and (jV 6 ) = 1.05. 
We see that after the transient regime (r > 15), which is 
proportional to the initial number of quanta, all curves 
present the same behavior - a linear time dependence 
with the same photons creation rate. In Fig. fc we 
plot the atomic population inversion (o~ z ) for these states, 
showing that asymptotically (a z ) approaches zero for any 
initial state. This result was expected due to the atomic 
decoherence. 

The master equation describes only the net effect of the 
environment on the system. For a better understanding 
of the role played by the atomic dephasing on the creation 
of photons we use the quantum trajectories approach |49[ 
I50j to study (n) during individual trajectories. Here the 
quantum jump operator is J p = "/ p h&zPO~ z and the non- 
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FIG. 1: Dynamics of the Rabi Hamiltonian for A = and 
g = 0.1 as function os dimensionless time r = ft (£ = 0.1). 
a) Mean photon number (n) for initial state \g, 0) without 
dephasing (line 1) and with dephasing 7 P h = 0.1 (line 2). 
There is photon generation from vacuum due to atomic de- 
phasing and the created field state is superpoissonian, since 
the Mandel factor q > (line 3). b) (n) for different initial 
states \4>i), i — 1, •■, 6 (see the text), demonstrating that the 
asymptotic photons generation does not depend on the initial 
state, c) Population inversion (a z ) for the states \<f>i) shown 
in (b): as expected, (<r z ) goes to zero asymptotically due to 
the decoherence. 



Hermitian Hamiltonian is H = H — i (j p h/2) I. In Fig. 
2 (a)-(c) we plot (n) versus r for 3 samples of individual 
trajectories for the initial number state Iff, 5), where we 
notice that in each trajectory (n) tends to increase as 
the time goes on. This occurs for two reasons: (i) the 
antirotating term in RH and (ii) the limitation of the 
cavity field from below by the vacuum, a |0) = 0. 

To illustrate (i) we plot in Fig. 2d (n) obtained via 
quantum trajectories approach for the JCH under the 
atomic dephasing, where we see that (n) oscillates with 
time but does not increase, contrary to the Figs. 2 (a)- 
(c). To explain the process of photon creation, we notice 
that any state may be written in terms of the basis states 
{|s, n)} (with s — {ff, e} and \n) the Fock state). Between 
the jumps, the system evolves according to the Rabi 
Hamiltonian ([!]) (the non-hermitian part is not important 
due to the normalization condition |51j ) The JC term 
[ff(a<7_|_ + a^<7_)] promotes \g, n) «-> \e,n — 1) transitions, 
while the antirotating term induces |ff, n) «-> |e, n + 1). 
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FIG. 2: a-c). Mean photon number for RH (n) versus t 
for three particular trajectories, obtained using the quantum 
trajectories approach. The parameters are A = 0, g = 0.1, 
7ph = 0.1 and the initial state is \g, 5). (n) tends to increase 
with time, d) (n) for a single trajectory using JCH, showing 
that there is no increase of the photon number without the 
antirotating terms, e) Asymptotic increase of mean photon 
number for the RH (line 1) and the test Hamiltonian He (line 
2, see text), demonstrating that the phenomenon is due to the 
antirotating terms and the limitation by vacuum, a|0) = 0. 



The combined action of both parts generates all the 
possible transitions. However, the Fock states are lim- 
ited from below by the vacuum state (ii), so there are 
more available states \s,m > n) than \s,m < n) and the 
mean number of photons between the jumps tends to in- 
crease. Upon a jump, the reservoir reads out the atom's 
state (through the application of a z on the wavefunction, 
which transforms |ff) — ► —^/j p h Iff) and |e) — > ^/J p h |e)) 
so the coherence between the states |ff) and |e) is lost and 
subsequent evolution under H will not bring the system 
to the state at the moment of the previous jump. For this 
reason after each jump (n) tends to be larger than upon 
the last jump, as clearly demonstrated in Figs. 2(a)-(c). 
After making a statistical average over many trajecto- 
ries one finds out that (n) always increases, in agreement 
with Fig. 1. 

One could suspect a third explanation for photons gen- 
eration due to atomic dephasing - the different weights 
y/n and y/n + 1 arising upon operating a and on the 
Fock state |n). To show that this is not the case we con- 
sider the test Hamiltonian He obtained through the sub- 
stitution of the operators a and in RH nip by the expo- 

— 1/2 

nential phase operators |52l 153] EL = (n + 1) a and 
EL, respectively, where E + E- = 1 — |0) (0|. By 
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FIG. 3: Asymptotic photons generation rate f3 = 
d (n(r)} /drlr^ao as function of a) (ui + luo) for fixed g e 7 P ^; 
b) g for fixed 7 P h and ojo; c) 7 P h for fixed g and t^o- We ob- 
serve that (3 ~ 7 P h and is inversely proportional to ujq +u). In 



the weak coupling regime (3 ~ c? 



doing this wc eliminate the weight factors and y/n + 1 
from the RH, since i? + |n) = |n + 1) and |rt) = 
(1 - S n0 ) \n - 1). In Fig. 2e we plot (n) for RH (line 1) 
and (he) obtained using He (line 2) versus r for initial 
state |0) = | (7, 0) . In both the cases there is photon cre- 
ation from vacuum, although (he) increases slower than 
(n). Therefore this photon creation phenomenon is not 
due to the different weights y/n and \Jn + 1 attributed 
to the operators a and , but to the presence of the EM 
vacuum state. 

From Fig. 1 we notice that asymptotically (n(r)) in- 
creases linearly with time, so one may gain a deeper in- 
sight into the problem by analyzing how the asymptotic 
photons generation rate (3 = d (n(r)) /d-rl^oo scales 
with u} + uj q , g, and 7^ (here we fix o> and vary other pa- 
rameters). In Fig. 3a we plot (3 versus [uj + ujq) for fixed 
g e jph, where we see that (3 increases when ujq decreases, 
so (3 is inversely proportional to (uj + ujo)- Fig 3b shows 
the dependence of (3 on g for fixed "/ p h and u>q : for g <C 1 , 
the weak coupling regime, the analysis of the curve shows 



that (3 ~ g 2 , however such a dependence is modified for 
large values of g. Finally, Fig. 3c shows that (3 ~ 7^ for 
fixed g and ujq, in agreement with the quantum trajec- 
tories approach: indeed, larger j p h implies larger jump 
probability and, consequently, in average (n) increases 
faster with time. A quantitative analysis of (3 will be 
presented elsewhere. 



III. RANDOM FREQUENCY FLUCTUATIONS 

One of the origins of dephasing from the physical point 
of view are the random shifts of the atomic transition fre- 
quency loq due to the interaction with the environment 
[4^1 IM] , Indeed, 1// noise in the bias controlling the 
atomic transition frequency is the dominant source of 
decoherence in superconducting artificial atoms (qubits) 
[3"2l [551 1551 157] . To investigate the effect of such a noise 
on the dynamics of the atom-cavity system, we integrated 
numerically the master equation for the RH (neglect- 
ing the damping and dephasing) assuming that uq has 
stochastic fluctuations. Our goal is to show that, when 
averaged over ensemble, this source of decoherence does 
asymptotically generate photons from vacuum, since its 
mathematical description is given by the dephasing Lind- 
blad superoperator in the master equation we studied 
above. 

As a simple model we considered time-dependent Wq(<) 

( O.lexr if cj (t) < fl - 0.8 e 
iv (t + dt) = u (t) + < -O.lsx r if u (t) > fl + 0.8 e , 
I O.lex (r — 1/2) otherwise 

(7) 

where n = if, = 0) is the mean atomic transition fre- 
quency, r € (—1,1) is a random number, e <C 1 is the 
maximum shift of the atom frequency and dt is the sim- 
ulation unit step, gdt <C 1. Here x is related to the 
'frequency' of the noise: qualitatively, for iClwe have 
'low frequency' noise, and in the opposite limit we have 
'high frequency' noise. 

We assumed the initial state \g, 0) and calculated the 
ensemble average of the mean photon number (n) av and 
the probability of exciting the atom P e using the param- 
eters Qq = lj 9 = 6 • 10~ 2 , e = g. We considered three 
examples of noise whose spectra are shown in Fig. 4a: 
x = 1 corresponds to the 'low frequency' noise (line 1) 
and x — 6 is our reference to the 'high frequency' noise 
(line 3), while x — 3 is a case in between we call 'mid- 
dle frequency' noise (line 2, data not shown in Fig. la 
since it lies between lines 1 and 3). First, the Fig. 4b 
shows three samples of (n) obtained for single runs of 
simulation for the 'high frequency' noise - one may see 
random rises and falls of photon number, however in av- 
erage (n) increases. The Figs. 4c and 4d show (n) and 
P e , respectively, obtained after averaging out many runs 
of simulations for the three kinds of noise. We see that 
in average there is a growth of both (n) and P e ; more- 
over, the growth of (n) av is approximately linear in time, 
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FIG. 4: Simulation of atomic dephasing via random atomic 
frequency fluctuations for parameters J7q = 1, g = 6 • 1CP 2 , 
e — g. a) Spectrum of the frequency noise u>o(t) — Qo. Line 
1 (red) corresponds to the 'low-' and line 3 (black) - to the 
'high-' frequency noises, b) (n) for three individual runs of 
simulation using high frequency noise, c) (n) averaged out 
over many runs of simulation, showing photons growth, de- 
pendent on the noise frequency. Here line 2 (blue) is the 
'middle-frequency' noise, d) P e averaged out over many sim- 
ulations. These curves agree qualitatively with the results 
obtained above using the master equation approach. 

in agreement with our previous results. 

From Fig. 4c we see that the photon generation rate 
is higher for higher frequency noise. One may under- 
stand qualitatively such a behavior as follows. The dy- 
namics of the RH with externally prescribed non-random 
ujo (t) allows for the coherent generation of both EM and 
atomic real excitations from vacuum for, e.g., periodic 
[HI BH S3 HE] or linear gTJ H3] time-dependence of 
wo (t). Moreover, the photon creation rate strongly de- 
pends on the shape of u>q (t) [4J] and, in the periodic case, 
on the periodicity of the modulation of u (t) HH US]- 
The Fourier transform (Fig. 4a) of the noise contains 
the 'resonant' frequencies (y ~ 2 j44j [46] ) , with respec- 
tive weights, for which photon generation occurs in the 
periodic case, so in average one expects a slow incoher- 
ent photon creation from vacuum due to these compo- 
nents in the noise spectrum. For the high frequency noise 
there are more resonant frequencies in the noise spec- 
trum and/or their weight is larger compared to the low- 
frequency noise, so the photon growth rate is higher, in 
agreement with Fig. 4c. 

Therefore, one of the physical origins of the photon 
creation through decoherence in cavity QED are the ran- 
dom fluctuations of the atomic transition frequency giv- 
ing rise to effective time-dependent RH, for which pho- 
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FIG. 5: a) (n) versus r using the master equation |2| for 
initial state \g,0) with parameters A = 0, g — 0.1, n t = 

and decay rates (7 P /i,7, k). Dotted line: (0,0,0), line 1: 
(0,0, 1) ■ 10" 1 , line 2: (0, 1,0) • HT 1 , line 3: (1, 1, 1) • 10"\ 
line 4: (1, 0, 1) ■ 10"\ line 5: (1, 1, 0) ■ 10" 1 . b) (n) versus gt 
for initial state \g,0) using circuit QED parameters A = 0, 
g = 2 ■ 10 -2 , nt = 6 • 10 -2 and distinct decay rates (line 

1 denotes the thermal photon number n t ). Line 2: current 
parameters (2,3,0.4) • 10~ 4 . Line 3: expected future scenario 
(2,3,0.4) • 10" 5 . Line 4: highly biased noise (200,3,0.4) • 
10 -4 . For current parameters there is no observable difference 
between nt and (n), however in the future or in very noisy 
environments (n) may become significantly larger than nt due 
to photon creation through decoherence. c) P e corresponding 
to the parameters in (b); P e resembles the behavior of (n) . 

tons are created from vacuum for the 'resonant' frequen- 
cies HH HB] present in the noise spectrum. 



IV. DEPHASING PLUS RELAXATION 

In realistic situations, besides the dephasing reservoir 
there are other important error sources (environments) 
acting on the system, e.g., thermal reservoirs. When 
other reservoirs are present, there is a competition be- 
tween photon creation due to the atomic dephasing and 
photon losses due to the damping. In order to see this 
effect, we plot in Fig. 5a (n) versus r for different decay 
rates in Eq. for the parameters A = 0, g = 0.1, 

n t = and the initial state \g,0). The effect of temper- 
ature (n t > 0) is just to shift the curves upward. The 
dotted line shows (n) in absence of any reservoir. When 
the atomic phase reservoir is switched off, the curve 1 
with dissipation parameters (0, 0, 1) • 10 _1 [we use nota- 
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tion (7 P /s,7, k)] and 2 with (0,1,0) • 10 _1 show that for 
large times (n) is smaller than in the cases where the 
phase reservoir is switched on, as shown in the curves 
3 with (1, 1, 1) • 10- 1 , 4 with (1,0, 1) • 10" 1 and 5 with 
(1, 1,0) • 10 _1 . Even in situations in which the environ- 
ment starts at T — OA", in the asymptotic limit the sys- 
tem behaves as being subjected to an effective reservoir 
with fit > 0, since lirriT-^oo (n(r)) > 0. The number of 
effective reservoir photons fit increases when the atomic 
phase reservoir is present (see curves 3, 4, and 5) and 
decreases when the atomic and field thermal reservoirs 
are predominant. However, the effective reservoir cannot 
be compared to the usual thermal reservoir, since the 
statistics of created field state is quite different from the 
thermal state statistics. A similar conclusion was drawn 
in |35l . where the authors considered the master equa- 
tion ([2]) at zero temperature with ^y p h = 7 = 0, k =/= 
and showed that the antirotating part of the RH gives 
rise to a thermal-like term in the effective master equa- 
tion, although it cannot be interpreted as an interaction 
with a thermal bath at a certain temperature. 

In fig. 5b we consider current experimental param- 
eters taken from recent circuit QED experiments [55] 
A = 0, g = 2 ■ 10~ 2 , n t = 6 • 10~ 2 . The line 1 shows 
the thermal photon number n t and the line 2 shows (n) 
obtained via master equation ^ for current dissipation 
rates: n = 4-10" 5 , 7 = 3-10~ 4 , j ph = 2-10~ 4 . There is no 
visible deviation of (n) from n t for the current tempera- 
tures, even if the detuning is increased to A = —0.2 (data 
not shown). Next we consider the scenario expected in 
the future, where the damping losses are suppressed one 
order of magnitude, k = 4 ■ 10 -6 , 7=3- 10~ 5 , but the 
dephasing rate remains the same, 7^ = 2 • 10~ 4 . In this 
case (line 3), (n) slightly deviates from the thermal pho- 
ton number n t and such an effect could be observed in 
very accurate measurements. Finally, we take the cur- 
rent values of damping rates, k = 4 • 10~ 5 , 7 = 3- 10~ 4 , 
and consider a highly biased noise |57j with the dephas- 
ing rate two order of magnitudes larger than the best one 
available today, 7^ = 2 ■ 10 -2 (line 4). In this case, (n) 
equals almost twice the thermal photons number due to 
the phenomenon of photon creation through decoherence. 

Finally, in Fig. 5c we plot P e corresponding to the 
parameters of Fig. 5b, where the line 1 shows p^ WA 
obtained for current dissipation parameters using RWA, 
which nearly does not depend on the relaxation rates. 
The behavior of P e resembles the one of (n) and indi- 
cates that P e increases due to the combined action of de- 
phasing and the antirotating term, although for current 
parameters this phenomenon is insignificant. However, 
for a large 7 P /j (line 4) P e is substantially higher than 
P^ WA , demonstrating that large dephasing, besides de- 
coherence, may also induce bit flip error. 



V. DISCUSSION AND CONCLUDING 
REMARKS 

We studied numerically the dynamics of the Rabi 
Hamiltonian subjected to dissipative losses, assuming ad 
hoc that one may describe the dissipative dynamics of the 
RH by using the standard master equation with usual 
damping and dephasing Lindblad superoperators. We 
found out that the atomic dephasing, when combined 
with the antirotating term in the RH, induces photon 
creation from vacuum. The physical interpretation of the 
phenomenon was given using two alternative approaches: 
1) the quantum trajectories approach based on quan- 
tum jumps and 2) microscopic ad hoc model of dephas- 
ing based on stochastic oscillations of the atomic transi- 
tion frequency (as occurs in solid state cavity QED). We 
showed that the photon creation through atomic deco- 
herence is suppressed in the presence of damping mecha- 
nisms, and estimated the magnitude of this phenomenon 
using currents experimental values of parameters, noting 
that the phenomenon might become relevant in future 
experiments. 

We do not have a formal proof that the master equa- 
tion ^ we used throughout this work is valid for the Rabi 
Hamiltonian in the strong atom-cavity coupling regime, 
since we assumed the Lindbladian dissipative superopera- 
tors without a microscopic deduction of the master equa- 
tion. Nevertheless, the results obtained here are relevant 
for two reasons. First, if the master equation ^ is in- 
deed applicable, it says that decoherence induces photons 
generation from vacuum, and such an effect may become 
relevant in future experiments with lower temperatures 
and lower damping rates. Second, if the master equation 
([2]) turns out to be non-applicable to this problem, it will 
call attention that there is an incompatibility between 
the standard Lindbladian superoperators and antirotat- 
ing terms in the Rabi Hamiltonians. In any case, more 
investigation of the problem is needed, since up to now 
the theoretical and numerical investigations were mainly 
concerned about the role of the RH antirotating term in 
the closed system dynamics, while our study points out 
novel important effects of the antirotating term in open 
systems. 
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